Cylindrical ideal magnetohydrodynamic equilibria 
with incompressible flows 

G. N. Throumoulopoulost and H. Tasso* 

t Section of Theoretical Physics, 
Physics Department, University of Ioannina, 
GR 451 10 Ioannina, Greece 

* Max-Planck-Institut fur Plasmaphysik, 
EURATOM Association, D-85748 Garching, Germany 



Abstract 

It is proved that (a) the solutions of the ideal magnetohydrodynamic equation, 
which describe the equlibrium states of a cylindrical plasma with purely poloidal 
flow and arbitrary cross sectional shape [G. N. Throumoulopoulos and G. Pantis, 
Plasma Phys. and Contr. Fusion 38, 1817 (1996)] are also valid for incompressible 
equlibrium flows with the axial velocity component being a free surface quantity 
and (b) for the case of isothermal incompressible equilibria the magnetic surfaces 
have necessarily circular cross section. 



I. Introduction 



In a recent paper |I[ it is proved that, if the ideal MHD stationary flows of a 
cylindrical plasma with arbitrary cross sectional shape are purely poloidal, they 
must be incompressible. This property simplifies considerably the equilibrium 
problem, i.e. it turns out that the equlibrium is governed by an elliptic partial 
differential equation for the poloidal magnetic flux function ip which is amenable 
to several classes of analytic solutions. For an arbitrary flow, i.e. when the ve- 
locity has non vanishing axial and poloidal components, the equilibrium becomes 
much more complicated. With the adoption of a specific equation of state, e. 
g. isentropic magnetic surfaces f2j, the symmetric equilibrium states in a two 
dimensional geometry are governed by a partial differential equation for ip, which 
contains five surface quantities (i.e. quantities solely dependent on ■?/>), in con- 
juction with a nonlinear algebraic Bernoulli equation. The derivation of analytic 
solutions of this set of equations is difficult. 

In the present note we study the equlibrium of a cylindrical plasma with 
incompressible flows and show that the incompressibility condition makes it pos- 
sible to construct analytic equilibria, which constitute a generalization of the ones 
obtained in Ref. |L|. This is the subject of Sec. II. The special class of incom- 
pressible equilibria with isothermal magnetic surfaces is examined in Sec. III. 
Section IV summarizes our conclusions. 

II. Equilibrium equations and analytic solutions 

The ideal MHD equilibrium states of plasma flows are governed by the fol- 
lowing set of equations, written in standard notations and convenient units: 



V ■ (pv) = (1) 

p(y. V)v = j x B- VP (2) 

V x E = (3) 
VxB=j (4) 

V • B = (5) 
E + vxB = 0. (6) 



The system under consideration is a cylindrical plasma with flow and arbitrary 
cross sectional shape. For this configuration convenient coordinates are £, r\ and 
z with unit basis vectors e^, e^, e z , where e z is parallel to the axis of symmetry 
and £, r\ are generalized coordinates pertaining to the poloidal cross section. The 
equilibrium quantities do not depend on z. The divergence free fields, i.e. the 
magnetic field B, the current density density j and the mass flow pv, can be 
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expressed in terms of the stream functions i/j(£,T]), F(£, 77), BziCiV) an d v z (£,r]) 
as 

B = B z e z + e z xVip (7) 
j = V 2 ^e z - e z x VB Z (8) 

and 

pv = pv z e z + e z x VF. (9) 

Constant ip surfaces are the magnetic surfaces. Eqs. ([[])- (0) can De reduced by 
means of certain integrals of the system, wich are shown to be surface quantities. 
To identify two of these quantities, the time independent electric field is expressed 
by E = — V$ and the Ohm's law (||) is projected along e z and B, respectively, 
yielding 

e z ■ (e 2 x VF) x (e z x Vip) = (10) 

and 

B-V$ = 0. (11) 

Eqs. (|10D and (0) imply that F = F{^) and $ = $(^). Two additional surface 
quantities are found from the component of Eq. (|6|) perpendicular to a magnetic 
surface: 

?a!--v z = &, (12) 
p 

and from the component of the momentum conservation equation (fj) along e z : 

B z -F'v z =X(^). (13) 

(The prime denotes differentiation with respect to ip). Solving the set of Eqs. 
( |12"1) and (|T3| ) for B z and v z , one obtains 
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and 



F'(ip)X(ip) - $'(^) 



(14) 



(15) 



p-(F'(V>)) 2 ^ 

With the aid of Eqs. (|T(1D-(|T3"D, the components of Eq. (|2|) along B and perpen 
dicular to a magnetic surface, respectively, are put in the form 



B 



VP 







(16) 
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and 



(F 



/\2 S 



+P 



IV# 



VPl 



0. 



(17) 



It is pointed out here that Eqs. (|l6|) and (|T^) are valid for any equation of state 
for the plasma. 

In order to reduce further the equilibrium equations, we employ the incom- 
pressibiliy condition 

V-v = 0. (18) 



Then Eq. (1) implies that the density is a surface quantity, 

P = pW, 

and, consequently, Eqs. ([14]) and ([15]) yield 

B z = B z (ip), v z = v z (ip). 



(19) 



(20) 

With the use of Eqs. ([H]) and (]20|) , Eq. fli~6|) can be integrated yielding an 
expression for the pressure, i.e. 



?/2 



2p 



(21) 



We note here that, unlike in static equilibria, in the presence of flow magnetic 
surfaces do not coincide with isobaric surfaces because Eq. (0) implies that B- VP 
in general differs from zero. In this respect, the term P s (ip) is the static part of 
the pressure which does not vanish when F' is set to zero; Eqs. fll4|) , ( |15"1) and 
(17) have a singularity when 



in 
p 



i. 



On the basis of Eq. @ for pv and the definitions 



(22) 



j Ap = — 7j — for the Alfven 
velocity associated with the poloidal magnetic field and the Mach number M 2 = 
4-, Eq. (f22|) can be written as M 2 = 1. 

v Ap 

Assuming now v ' ^ 1, and inserting Eq. (|21[) into Eq. (]17|), the latter 
reduces to the elliptic differential equation 



1 - 



(F 



r\2 



P 



v 2 v + 



7 \2~P 



-F" |V^| 2 + [P s + 



Bi 



(23) 
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The absence of any hyperbolic regime in Eq. fl23|) can be understood by noting 
that, as is well known from the gas dynamics, the flow must be compressible to 
allow the equilibrium differential equation to depart from ellipticity. Eq. ( |23| ) 
does not contain the axial velocity v z and is identical to the equation governing 
cylindrical equilibria with purely poloidal flow jjfl. With the use of the ansatz 
d F" (F 1 ) 2 

t-p = 2-p-, which implies that v ^ = M 2 = const., Eq. Q2BD reduces to 

v2 "' + T^W( P ' + f)' = a (24) 

This is similar in form to the equation governing static equilibria; the only explicit 
reminiscence of flow is the presence of M c . Eq. fl24|) can be linearized for several 
B 2 

choices of P s + and a variety of analytic solutions of the linearized equation 
can be derived. In particular, the exact solutions for a circular cylindrical plasma 
obtained in Ref. ]IJ are also valid for incompressible equilibrium flows with a free 
axial velocity v z (ip). 

The singularity M 2 = 1 is the limit at which the confinement can be assured 
by the axial current V 2 ip alone. For M 2 > 1 the derivative of B 2 /2 must partly 
compensate for the pressure gradient. 

III. Equilibria with isothermal magnetic surfaces 

For fusion plasmas the thermal conduction along B is fast compared to the 
heat transport perpendicular to a magnetic surface and therefore equilibria with 
isothermal magnetic surfaces are of particular interest. The plasma is also as- 
sumed to obey the ideal gas law P = RpT. For this kind of equilibria, Eq. (|2T| ) 
implies that |V^| is a surface quantity and consequently from Eq. (fi7|) it turns 
out that V 2 ip is a surface quantity as well. Thus, the incompressible, T = T(^) 
equilibria satisfy the set of equations 

I W| 2 = {gm 2 (25) 

and 

V 2 ^ = /(# (26) 

Eqs. (^) and (|26|) imply that, on a magnetic surface the modulus of the vector 
V"0, which is perpendicular to this (arbitrary) magnetic surface, and V 2 ip, related 
to the variation of \Vip\, are constants. Therefore, one could speculate that 
magnetic surfaces are restricted to be circular. This conjecture can be proved as 
follows. 
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The coordinates £, rj and z are specified to be the Cartesian coordinates x, y, 
z. With the introduction of the quantities p = dip/dx, q = dip/dy, r = d 2 ip/dx 2 
and t = d 2 ip / dy 2 , Eqs. fl25|) and (EH) are written in the form 



and 



2,2 2 

P +Q =9 



r + t = f. 



(27) 



(28) 



The set of Eqs. (p7|) and ( ]28|) can be integrated by applying a procedure sug- 
gested by Palumbo || . Accordingly, considering the functions p and g which are 
functions of x and ?/ as functions of x and ip(x, y) one has 



<9p 

dx 



dp dp 
dx ^ dip 



and 



t ^ 
^ dip' 



(29) 



(30) 



(It is noted here that a surface function ( = ((x,y) = C(V ; ) can De employed 
instead of With the aid of Eqs. (|2?|), (|2|) and (gUp, Eq. @ reduces to 

<5x , 



/ — gg' and consequently 



P = x(f-gg')+h(iP). (31) 
On a magnetic surface it holds that dip = j^dx + ^jdy = 0, and therefore 



dy 
dx 



p_ = 
g 2 ^2 



[x (f - gg') + hf 



[x(f-gg') + hy 



(32) 



Introducing the new quantities a(ip) = f — gg', X = ax + h and Y = ay, Eq. 
(|32"D is put in the form 



dX 



X s 



(33) 



g 2 -X 2 ' 

Eq. (|33| ) describes a circle on the (x, y) plane with radius \g\ centred at (—h/a, 0) 

IV. Conclusions 



It was proved that the ideal MHD equilibrium states of a cylindrical plasma 
with incompressible flows and arbitrary cross section shape satisfy an elliptic 
partial differential equation [Eq. fl2"3"|)1, which is identical to the equation gov- 
erning cylindrical equilibria with purely poloidal flow; the axial flow velocity is 
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a free surface quantity. This equation permits the construction of several classes 
of analytic solutions. In particular, the exact equlibrium solutions for a circular 
cylindrical plasma and purely poloidal flow |T[ are also valid for the present case. 
In addition, it was proved that the magnetic surfaces of isothermal incompressible 
equilibria must have circular cross section. 

It is intersting to investigate symmetric incompressible equlibria in geometries 
representing more realistically the magnetic confinement systems, e.g. axisym- 
metric and straigth helically symmetric configurations. In this respect it may 
be noted here that, as proved in Ref. [f§], the special class of axially symmetric, 
incompressible, /3 P = 1, MHD equilibria with purely poloidal velocity does not 
exist; the only possible stationary equilibria of this kind are of cylindrical shape. 
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